OPEN Citation: Oncogenesis (2014) 3, e86; doi:1 0.1 038/oncsis.201 3.48 flp*. 

© 2014 Macmillan Publishers Limited All rights reserved 2157-9024/14 



www.nature.com/oncsis 



ORIGINAL ARTICLE 

Tumor suppressor NDRG2 tips the balance of oncogenic 
TGF-p via EMT inhibition in colorectal cancer 

L Shen 1 ' 7 , X Qu 2 ' 7 , Y Ma 1 ' 7 , J Zheng 3 , D Chu 1 , B Liu 4 , X Li 1 , M Wang 5 , C Xu 3 , N Liu 6 , L Yao 1 and J Zhang 1 

Transforming growth factor-beta (TGF-P), a pluripotent cytokine expressed in the colon, has a crucial but paradoxical role in 
colorectal cancer (CRC). TGF-p is a potent proliferation inhibitor of normal colon epithelial cells and acts as a tumor suppressor. 
However, TGF-p also promotes invasion and metastasis during late-stage CRC, thereby acting as an oncogene. Thus, understanding 
the factors behind the paradoxical roles of TGF-p and elucidating the mechanisms by which TGF-p-induced proliferation inhibition 
is impaired in CRC are necessary. Here, we found that the N-Myc tumor suppressor gene downstream-regulated gene NDRG2 
(N-Myc downstream-regulated gene 2), which is a TGF-p-responsive gene, abrogated TGF-p-induced epithelial-mesenchymal 
transition (EMT) and further inhibited the invasion and migration of CRC cells. TGF-p positively induced NDRG2 expression through 
direct transactivation mediated by Spl and by abrogation of the repressive c-Myc/Miz-1 complex on NDRC2 promoter in normal 
epithelial cells. Aberrant hypermethylation of NDRC2, which could respond to TGF-p growth inhibition signaling, abrogated the 
inhibitory effect of NDRG2 in TGF-p-induced EMT in CRCs. Reduced NDRG2 expression was highly correlated with the invasion stage 
and metastasis of CRC. Our study establishes that NDRG2 is a new tumor suppressor gene that responds to TGF-p anti-proliferative 
signaling and tips the balance of oncogenic TGF-p during late-stage CRC. 
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INTRODUCTION 

Colorectal cancer (CRC) is a commonly diagnosed cancer in men 
and women, and it is the third most common fatal malignancy 
and second leading cause of cancer-related deaths in the western 
world. 1 The etiology of CRC is multifactorial, and it is widely 
accepted that genetic factors, such as Ras/Raf, transforming 
growth factor (TGF)-p, p-catenin/T-cell factor and the p53 
pathway, have key roles in the predisposition to cancer 
development and progression. 2 Although genetic changes have 
been identified at multiple levels throughout the multistage 
carcinogenesis process, the initial factors corresponding to these 
changes need to be further elucidated. 

TGF-p is a secreted multifunctional cytokine that regulates a 
variety of cellular processes. The TGF-p receptor complex activates 
the Smad signaling pathway, which includes Smad2, Smad3 
and Smad4, and non-Smad signaling cascades (for example, PI3K/ 
AKT, p38MAPK, MAPK-ERK and JNK) to produce the full spectrum 
of TGF-p responses. 3 5 The role of TGF-p signaling in 
carcinogenesis is complicated. Originally named for its 
transforming activity in in vitro assays, TGF-p now unequivocally 
demonstrates both tumor suppressor and oncogenic activities. 
The tumor suppressor activities dominate in normal tissue and 
mainly occur through the direct regulation of cell-cycle inhibitors, 



such as p21 cipl and p i5 INK4B , 6 ' 7 and cell-cycle activator c-Myc via 
transcriptional and post-transcriptional mechanisms. 8 However, 
during tumorigenesis, changes in TGF-p expression and cellular 
responses tip the balance in favor of oncogenic activities by 
inducing the epithelial-mesenchymal transition (EMT), which is 
mediated by Fibronectin, Twist, Snail and so on, and finally 
accelerating tumor invasion and metastasis. 911 There is 
considerable genetic evidence that the loss of sensitivity to 
growth inhibition by TGF-p is an important event in colorectal 
carcinogenesis. Much of the evidence is derived from studies in 
human CRCs demonstrating inactivating mutations in genes 
encoding proteins involved in TGF-p signal transduction, 
including SMAD4, 12 SMAD2, 13 and TGFBR2. 14 However, it has 
also been reported that restoration of an impaired TGF-p pathway 
cannot restore the anti-proliferative response to TGF-p in CRC 
cells. 15,16 Therefore, to fully understand the paradoxical effect of 
TGF-p in carcinogenesis, other factors and mechanisms need to be 
uncovered and elucidated. 

In recent years, a new tumor suppressor gene, NDRG2 (N-Myc 
downstream-regulated gene 2), belonging to the NDRG family that 
consists of four identified members (A/DRG1-4), has been 
implicated in the regulation of cell differentiation and prolifera- 
tion. The role of NDRG2 in cancer was supported by the initial 
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observation that NDRG2 inhibited the proliferation of glioblastoma 
cells 17 and promoted the differentiation of various cell types, 
including dendritic cells, 18 PCI 2 cells 19 and CRC cells. 20 Consistent 
with the potential function of NDRC2 as a tumor suppressor, 
downregulation of NDRC2 expression was found in multiple 
human cancer cell lines and tumors, particularly colorectal 
carcinomas and high-risk adenoma precursor lesions. 21 Most 
notably, our recent data showed that a decrease in NDRG2 
expression in primary human CRC is a powerful and independent 
predictor for poor prognosis of CRC, particularly those with 
advanced tumor node metastasis (TNM) stage and lymph-node 
metastasis. 22 Thus, the relevance of NDRG2 to the paradoxical 
roles of TGF-p in CRC should be investigated. 

In addition, accumulating evidence indicates that in addition to 
genetic changes, epigenetic alterations have a major role in the 
initiation and progression of CRC. 23,24 More importantly, tumor 
suppressor genes inactivated by mutations and promoter CpG 
island methylation are more frequently affected by 
hypermethylation than mutations in CRC. 25 At the single-cell 
level, mutation/copy number changes and hypermethylation are 
most frequently mutually exclusive 25 In addition, at the genome 
level, there is an inverse relationship between frequent epigenetic 
changes, the CpG island methylator phenotype, large-scale 
genomic gains or losses and chromosomal instability, 26,27 
which has been characterized in ~85% of CRCs with increased 
chromosomal losses and heterozygosity 28 30 

Our study demonstrates for the first time that NDRC2 is a TGF-p 
responsive gene that abrogates TGF-p-induced EMT in CRC 
cells. TGF-p upregulated NDRC2 via transcriptional activation 
mediated by the transcription factor Spl and by abrogation of the 
repressive c-Myc/Miz-1 complex in epithelial cells. Furthermore, 
aberrant NDRG2 hypermethylation, which could respond to TGF-p 
growth inhibition signaling, abrogated the inhibitory effect of 
NDRG2 in TGF-p-induced EMT in CRCs. Most importantly, the 
reduced NDRG2 expression was highly correlated with the 
invasion stage and metastasis of CRC. Together, our data 
show that NDRG2, as a TGF-p responsive tumor suppressor gene, 
tips the balance of the TGF-p role toward oncogenesis in 
late-stage CRCs. 



RESULTS 

Positive correlation between NDRG2 expression and the Smad2/3 
phosphorylation level in CRC samples 

Our group recently demonstrated that NDRG2 is an important 
tumor suppressor that correlates with metastasis-free and overall 
survival in CRC patients. 22 Because accumulated data suggest that 
the phosphorylation of Smad2/3 at its C-terminal region correlates 
with TGF-p-induced cell-cycle arrest and growth inhibition, 31,32 
we analyzed phosphorylated Smad2/3 levels to evaluate the 
tumor suppressive activity of TGF-p. Immunohistochemical (IHC) 
analysis of 50 pairs of human colorectal tumors and adjacent 
normal colon samples revealed that NDRG2 and phosphorylated 
Smad2/3 were highly expressed in the adjacent normal colon 
tissues compared with the colorectal tumor samples and were 
inversely correlated with tumor stage (Figure la). Moreover, 
phosphorylated Smad2/3 was predominantly located in nuclei, 
suggesting a role of Smad complex in transcriptional regulation. 
Indeed, in normal versus cancer tissue, p-Smad2/3 levels were 5.1 1 
versus 0.65 (Figure 1 b), respectively, and NDRG2 levels were 5.79 
versus 0.65 (Figure 1c), respectively. Additionally, even less 
staining was observed in the more aggressive CRCs (Tables 1 
and 2). Furthermore, a correlation study determined that NDRG2 
expression was positively correlated with phosphorylated Smad2/ 
3 expression in these tissue samples (P<0.001, ANOVA test) 
(Figure Id), which suggests that NDRG2 might be involved in 
TGF-p downstream signaling. 



NDRG2 is a new direct responsive target of TGF-p 
To further determine the regulation of NDRG2 by TGF-p, we first 
confirmed the inhibition of cellular proliferation and cell-cycle 
arrest due to TGF-p in the human keratinocyte cell line HaCaT, 
which has an intact TGF-p-Smad signaling pathway (Supplementary 
Figures SI A and B). Intriguingly, the NDRG2 protein and messenger 
RNA (mRNA) expression levels were increased in the presence of 
TGF-p (Figures 2a and b), which indicates the transcriptional 
regulation of NDRG2 by TGF-p. Since it has been shown that c-Myc 
has a key role in NDRG2 inhibition 20 and TGF-p-induced cell-cycle 
arrest, 33 we observed a decrease in the expression of c-Myc and 
increased NDRG2 expression in the presence of TGF-p (Figures 2a 
and c). In addition, TGF-p-induced growth inhibition was plotted 
against NDRG2 and the control shRNA (Figure 2d). A key step in 
defining direct transcriptional regulation is determining the core 
region within the NDRG2 promoter responsible for TGF-p induction. 
Thus, a series of previously created processive NDRG2 promoter 
truncation constructs 20 were used in the present study. The 
presence of TGF-p increased the full-length NDRG2 promoter 
construct (- 1455/+274) activity >3-fold, and even a 1.6-fold 
induction was observed when using a minimal core promoter 
construct (—79/ + 57) (Figure 2e). Finally, a chromatin immuno- 
precipitation (ChIP) assay was performed to validate the direct 
binding of the TGF-p downstream targets Smad2/3 and Smad4 
with the NDRG2 promoter, which was enhanced by TGF-p 
(Figure 2f). In contrast, mutation of the Smad binding site showed 
little effect on NDRG2 promoter activity (data not shown), 
indicating that other transcription factors might be involved in 
TGF-p-mediated NDRG2 induction. Thus, the NDRG2 tumor 
suppressor gene was determined to be a new downstream target 
gene in the TGF-p pathway in epithelial cells that responds to TGF-p 
growth inhibitory effects. 

Spl is necessary for the transcriptional activation of NDRG2 by 
TGF-p 

To fully understand the regulatory mechanism of NDRG2 by TGF-p, 
we analyzed the structure of the NDRG2 promoter, which was 
considered necessary and sufficient for induction by TGF-p. Three 
potential Spl binding sites were identified that are GC-rich motif 
in NDRG2 (-250/ +274) (Supplementary Figure S2). Intriguingly, 
Spl could upregulate NDRG2 at the protein and mRNA levels 
(Figure 3a). NDRG2 promoter activity was also induced by Spl 
overexpression (Figures 3a and b). However, Spl induction was 
dramatically reduced on NDRG2 promoter mutants that had the 
first two Spl consensus sites deleted or mutated (Figure 3b), 
suggesting that these two sites are critical for the Spl regulation 
of NDRG2 promoter activity. This binding activity was also 
confirmed by EMSA (Supplementary Figure S3) and ChIP assays, 
which showed the interaction between Spl and the NDRG2 core 
promoter region ( — 215 to —53 and —75 to +50) (Figure 3c). 
Moreover, the presence of TGF-p and overexpression of Spl 
synergistically induced a greater than 10-fold increase in the 
promoter activity of wild-type NDRG2 ( - 1 455/ + 274 and - 273/ 
+ 274), and there was only a slight influence on the mutant 
NDRG2 (- 1455/ + 274) promoter (Figure 3d). Accordingly, an 
endogenous Spl knockdown via shRNA (psilencer-Spl) or 
functional blocking using a dominant-negative Spl (pEBG-DN- 
Spl) decreased NDRG2 expression (Supplementary Figures S4A-C) 
and abolished TGF-p responsiveness either on wild-type NDRG2 or 
on Spl consensus sites mutant constructs in advance (Figure 3e). 
In addition, Mithramycin A, an inhibitor that inhibited Spl binding 
with DNA by modifying GC-rich sites, dose dependently reduced 
NDRG2 expression and promoter activity (Supplementary 
Figures S4D and E) and blocked TGF-p-induced NDRG2 activity 
(Figure 3f). Therefore, Spl binding to consensus sites in the NDRG2 
promoter was required for the TGF-p-mediated transcriptional 
activation of NDRG2. 
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Figure 1. Positive correlation between NDRG2 and p-Smad2/3 in colorectal cancer tissues, (a) IHC staining of NDRG2 and p-Smad2/3 in human 
adjacent normal and cancerous colon samples from patients at different stages. Case 1 is from a stage III patient, and Case 2 is from a stage II 
patient. Original magnification: x 5, x 40 or x 100. (b, c) Relative expression level of p-Smad2/3 and NDRG2 in human tissue samples. 
Student's f-test was applied for statistical analyses. *P<0.05, which was considered significantly different, (d) Correlation between NDRG2 and 
p-Smad2/3 expression with linear regression and Pearson's correlation significance (P< 0.001, ANOVA test). 



Table 1. Statistical results of the immunohistology 




n 


C-terminal p-Smad2/3 
+ to + + + 


P-value 


Total 


50 


34 


16 




Sex 








0.609' 


Male 


26 


18 


8 




Female 


24 


16 


8 




Age at diagnosis, years 








0.805' 


<60 


25 


16 


9 




>60 


25 


18 


7 




Differentiation status 








<0.001 b 


Poor 


9 


7 


2 




Moderate 


22 


17 


5 




Well 


19 


10 


9 




TNM stage 








<0.001 b 


l-ll 


8 


1 


7 




III 


16 


12 


4 




IV 


26 


21 


5 




Abbreviation: TNM, tumor node metastasis, 
levels were compared using Fisher's exact test, 
levels were compared using Kruskal Wallis. 


'P-value when expression 
b P-value when expression 



Table 2. Statistical results of immunohistology 




n 




NDRG2 


P-value 








+ to + + + 


Total 


50 


38 


12 




Sex 
Male 
Female 


26 
24 


20 
16 


6 
8 


0.147' 


Age at diagnosis, years 
<60 
>60 


25 
25 


18 
18 


7 
7 


0.173' 


Differentiation status 
Poor 

Moderate 
Well 


9 
22 
19 


8 
18 
10 


1 

4 

9 


<0.001 b 


TNM stage 
l-ll 


8 


1 


7 
5 
2 


<0.00l b 


III 
IV 


16 
26 


11 
24 




Abbreviations: NDRG2, N-Myc downstream-regulated gene 2; TNM, tumor 
node metastasis. a P-value when expression levels were compared using 
Fisher's exact test. b P-value when expression levels were compared using 
Kruskal Wallis. 



TGF-)3 abrogates the c-Myc/Miz-1 repressive complex to activate 
NDRG2 

On the basis of the above data, we observed that the removal of 
Spl consensus sites could not completely abolish TGF-p activity 
on the NDRC2 promoter (Figures 3d and e), suggesting that 



other factors contributed to the induction of NDRC2 by TGF-p. 
Our previous data revealed that c-Myc could transcriptionally 
repress NDRG2 by interacting with Miz-1 on the NDRG2 
promoter. 20 We next investigated the importance of the c-Myc/ 
Miz-1 complex in response of NDRG2 to TGF-p. Because c-Myc is a 
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Figure 2. TGF-p directly transactivates NDRG2 at the transcriptional level, (a, b) HaCaT cells were treated with TGF-p for 24 h, and the protein 
and mRNA expression levels of NDRG2 and p21 Qpl were detected. p-Actin served as a control to ensure equal loading, (c) HaCaT cells were 
treated with TGF-p for the indicated time, and the c-Myc and NDRG2 protein levels were detected. p-Actin served as a control to ensure equal 
loading, (d) HaCaT cells were treated with 10 ng/ml TGF-p for 48 h. Relative cell viability was detected with an MTT assay {n = 3). (e) A panel of 
NDRG2 promoter constructs were transfected into HaCaT cells. At 24 h post transfection, the cells were treated with 10 ng/ml TGF-p for 20 h. 
The data shown are the means ± s.d. from triplicate analysis, (f ) Top: Representation of the human NDRG2 promoter regions depicting the 
sequence of the Smad-binding region (SBR). Bottom: HaCaT cells after TGF-p treatment were extracted and subjected to a chromatin 
immunoprecipitation assay. The presence of the NDRG2 promoter was detected with RT-PCR and visualized with ethidium bromide staining. 



key transcription node during EMT 34 that has an important role in 
the TGF-p-induced expression of cyclin-dependent kinase 
inhibitors, we observed that in addition to the induction of 
NDRG2, TGF-p caused a rapid reduction in c-Myc (Figures 2a and 
c). Moreover, c-Myc overexpression dramatically abolished the 
NDRG2 and p21 Clpl expression induced by TGF-p (Figure 4a) and 
completely abolished TGF-|3-mediated induction of the NDRG2 
promoter, which contains mutant Spl binding sites (Figure 4b). In 
addition, there was TGF-p-induced NDRG2 promoter activity in 
cells overexpressing Miz-1, and this effect was further enhanced 
with overexpression of Miz-1 A(637-803), which lacks the C-term- 
inal residues (637-803) responsible for interaction with c-Myc 35 
(Figure 4c). However, c-Myc restoration counteracted the Miz-1 
activating effects of the NDRG2 promoter (Figure 4d) but failed to 
block Miz-1 A(637-803) activity (Supplementary Figure S5). 
Furthermore, treatment of cells with TGF-p did not dramatically 
alter the DNA binding activity of Spl on either of the Spl 
consensus sites (Figure 4e), whereas Smad4 and Smad2/3 binding 
to the Spl consensus sites was increased, indicating that Sp1 
could mediate Smad complex binding on the NDRG2 promoter. In 
addition, TGF-p stimulation significantly reduced the binding of 
c-Myc at the NDRG2 promoter but did not alter Miz-1 occupancy 
(Figure 4f). Collectively, our results demonstrated that TGF-p 
upregulated NDRG2, and that is occurred partially through the 
abrogation of the interaction between c-Myc and Miz-1. 

NDRG2 overexpression attenuates TGF-p-induced EMT 
In agreement with previous reports, 36,37 we observed that CRC 
cell lines with multiple TGF-p/Smad pathway mutations lose 
TGF-p-induced growth inhibition characteristics (Supplementary 
Figure S6). To further determine whether NDRG2 is involved in the 
TGF-p-mediated induction of EMT in CRC cells, we established 



stable HT29 cell lines overexpressing NDRG2 or p21 Clp1 . There 
were no obvious morphological differences in HT29 cells after 
NDRG2 overexpression, which is likely because the parental cells 
already exhibit an epithelial phenotype. After TGF-p treatment, 
HT29 cells with NDRG2 overexpression showed an epithelial 
morphology (that is, the cells clustered together in groups, 
compared with the fibroblast-like shape alteration in the control 
group) (Figure 5a). As determined by quantification, the reduction 
in the circularity indices was significantly reversed by NDRG2 
overexpression in the presence of TGF-p (Figure 5b). In addition, a 
reduction in E-cadherin and an increase in vimentin were observed 
in TGF-p-treated control cells. In contrast, NDRG2 overexpression 
inhibited the E-cadherin and vimentin alterations induced by TGF- 
P (Figures 5c and d). Moreover, NDRG2 overexpression decreased 
the migratory and invasive potential of HT29 cells and further 
reduced TGF-p-stimulated migration and invasion (Figures 5e and 
f). Notably, in HT29 cells overexpressing p21 Clp1 , which is one of 
the classical TGF-p responsive targets, EMT inhibition was not 
significant (Figures 5a-f). These studies indicated that NDRG2 could 
abrogate the migration, invasion and tumor progression induced 
by TGF-p in late-stage CRC via EMT inhibition. 



NDRG2 knockdown promotes EMT by upregulating the key factors 
involved in TGF-p-induced EMT 

To further confirm the EMT inhibitory ability of NDRG2, we 
established stable HT29 cell lines with either NDRG2 or p21 cip1 
knockdown and analyzed cell motility and invasion. The cell 
morphology of HT29 cells with the NDRG2 knockdown converted 
to a more elongated, fibroblast-like shape (Figures 6a and b), and 
TGF-p stimulation enhanced this morphological change. In 
addition, NDRG2 knockdown significantly induced EMT marker 
alterations (that is, vimentin increase and E-cadherin loss), which 
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Figure 3. TGF-p-induced NDRG2 induction requires the Spl transcription factor, (a) Cells were transfected with pcDNA3 or pcDNA3-HA-Sp1 . 
At 12 or 24 h after transfection, protein was extracted for western blotting (top), quantitative real-time PCR (middle) and a reporter assay 
using NDRG2 (— 1455/ + 274) (bottom), (b) HaCaT cells were transfected with the indicated NDRG2 reporter constructs with or without a 
Spl expression vector. Luciferase activities were detected 24 h after transfection. (c) Human intestinal epithelial cells (HIEC) transfected with 
pcDNA3 or pcDNA3-HA-Sp1 were immunoprecipitated with the indicated antibodies. The presence of the NDRG2 promoter containing Spl 
binding sites was detected by RT-PCR and visualized with ethidium bromide staining, (d, e) HaCaT cells were transfected with the indicated 
NDRG2 reporter constructs with or without Spl expression or shRNA vector. The cells were incubated in the presence or absence of 10 ng/ml 
TGF-(3 for 20 h before lysis and then analyzed for luciferase activity, (f) Quantitative real-time PCR (left) and luciferase reporter assay (right) 
analysis for the involvement of mithramycin in TGF-p-mediated NDRG2 expression. HaCaT cells were pretreated with or without mithramycin 
for 1 h at the indicated concentrations and then stimulated with 10 ng/ml TGF-p for 20 h. (a-f) The data are presented as the means ±s.d. 
from triplicate analyses. 



became more pronounced after the addition of TGF-p (Figure 6c). 
As expected, NDRG2 knockdown also significantly increased the 
migration and invasion of HT29 cells, and the potential for these 
processes was enhanced after incubation with TGF-p (Figures 6d 
and e). Similar results were observed in p21 Clpl knockdown cells 
(Figures 6a-e). Moreover, by screening for key factors involved in 
TGF-p-induced EMT, we observed that NDRG2 knockdown 
dramatically upregulated the fibronectin, Twist, Snail and N-cad- 
herin expression levels (Figure 6f). Together, these studies 



demonstrate that NDRG2 knockdown could promote EMT altera- 
tions in HT29 cells via the upregulation of multiple factors 
essential for TGF-p-induced EMT. 

Reduced NDRG2 expression caused by promoter 
hypermethylation fostered the progression and metastasis of CRC 

NDRG2 exhibited dramatically reduced expression at the mRNA 
and protein levels in the majority of CRC cells relative to normal 
human intestinal epithelial cells (HIEC) (Figure 7a). To determine 
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in the presence or absence of lOng/ml TGF-p, and the mRNA levels and luciferase activity were determined. The data are presented 
as the means ± s.d. from triplicate analyses, (a) Cells were transfected with either a control or a c-Myc expression vector. Quantitative real-time 
PCR was performed, (b) Cells were transfected with the indicated NDRG2 reporter constructs with or without the c-Myc expression vector. 

(c) Cells were transfected with the NDRG2 reporter construct together with wild-type Miz-1 and its mutant construct, as indicated. 

(d) Cells were transfected with the NDRG2 reporter construct with or without a c-Myc expression plasmid. Increasing amounts (0.01, 0.05 
and 0.1 |ig) of Miz-1 construct were co-transfected as indicated, (e, f) HaCaT cells treated with TGF-p were extracted and subjected 
to chromatin immunoprecipitation. The presence of the NDRG2 promoter was detected by RT-PCR and visualization by ethidium 
bromide staining. 



the mechanism for the resistance to growth inhibition and 
inability of TGF-p to induce NDRG2 in CRC, we performed DNA 
methylation analysis of the NDRG2 promoter in CRC cells and 
patient samples. Hypermethylation of the NDRG2 promoter was 
detected in most CRC cell lines (Figure 7b). Moreover, ~75% (18/ 
24) of the clinical colorectal tumors showed hypermethylation at 
the NDRG2 promoter, and a low percentage of hypermethylation 
(6/24) was observed in the normal counterparts (Figure 7c). 

To determine whether promoter hypermethylation in CRC leads 
to the reduced sensitivity of NDRG2 to multiple growth inhibition 
signals, including TGF-p, we analyzed the NDRG2 expression 
induced by Spl in NDRG2 hypermethylated CRC cells. Intriguingly, 
NDRG2 induction by Spl was completely abolished in NDRG2 
hypermethylated CRC cells (Supplementary Figure S7). Further- 
more, Spl binding activity was blocked by hypermethylation of 
the CpG dinucleotide within its binding site (A2) or adjacent CpG 
sites (A4) (Figure 7d). As expected, the DNA methyltransferase 
inhibitor 5-aza-2'-deoxycytidine (5-aza-dC) restored NDRG2 
expression, promoted the induction of NDRG2 by Sp1 
(Figure 7e), and further inhibited CRC proliferation 
(Supplementary Figure S8), thus indicating that epigenetic 
alteration is a critical step during tumor progression. However, 
the inhibition of EMT by 5-aza-dC was not dramatic (data not 
shown), which might be attributed to the global effects of 5-aza- 
dC on multiple target genes. Collectively, NDRG2 promoter 
methylation abrogated the Spl -meditated induction of NDRG2, 
which may be responsible for the reduction in the NDRG2 
response to TGF-p in CRC. 

On the basis of the above data, the NDRG2 tumor suppressor 
has an important role in impairing the tumor progression induced 



by TGF-p in CRC cells through inhibition of EMT. To further 
understand the role of A/D/?G2-mediated EMT inhibition in 
patients, we analyzed NDRG2 expression status and tumor 
invasion and metastasis characteristics in 260 CRC patients 
(Table 3). Consistent with the in vitro data, CRC patients with 
reduced NDRG2 expression showed a relatively higher invasion 
stage than those with a normal NDRG2 expression level (Figure 7f). 
Most importantly, reduced NDRG2 was significantly correlated 
with lymph-node involvement and distant metastasis (Figures 7g 
and h), confirming that NDRG2 is a potent tumor suppressor that 
acts not only through anti-proliferation effects but also via the 
inhibition of invasion and metastasis. 

As described in the working model (Figure 7i), our study 
identified NDRG2 as a new tumor suppressor gene that is a TGF-p- 
responsive target that tips the balance of oncogenic TGF-p via 
EMT inhibition in late-stage CRCs. 

DISCUSSION 

TGF-p serves as a tumor suppressor primarily by inhibiting cell 
proliferation during early tumor stages, and it behaves as a tumor 
promoter by stimulating both invasion and metastasis during late 
tumor stages. There is growing evidence that restoration of the 
impaired TGF-p pathway cannot rescue the anti-proliferative 
response to TGF-p in CRC cells. 9,16 Therefore, other mechanisms 
must be elucidated to determine how selective alterations in the 
TGF-p signaling pathway contribute to the development, 
progression and metastasis of CRC. 

In this study, we first demonstrated a clinical correlation 
between the protein accumulation of NDRG2 and C-terminal 
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Figure 5. NDRG2 overexpression attenuates TGF-p-induced EMT. (a) Phase-contrast photomicrographs of control cells and cells treated with 
lOng/ml TGF-p for 48 h. (b) Circularity index of untreated and TGF-p-treated cells, (c) Western blot analysis of E-cadherin and vimentin 
expression from total lysates of untreated cells and cells treated with 10ng/ml TGF-p for 48 h. (d) Relative mRNA levels of E-cadherin and 
vimentin in control cells and cells treated with 10 ng/ml TGF-p for 48 h. The control values were normalized to 1, and the data are expressed as 
fold-change in treated cells, (e, f) The migratory and invasive behavior of untreated cells and cells treated with 10 ng/ml TGF-p for 48 h. 
The control values were normalized to 1, and the data are expressed as fold-change in treated cells, (d-f) The data are presented as the 
means ±s.d. from triplicate analyses. 



phosphorylated Smad2/3, which responds to TGF-p anti- 
proliferative signaling, by the IHC staining of colon cancers; these 
two proteins are also negatively correlated with colon 
cancer progression. Our analysis of 50 pairs of colon cancer and 
adjacent normal colon samples demonstrated that TGF-p- 
mediated activation likely serves as an important guardian for 
the expression of NDRC2 in vivo. It should be noted that, as is 
often observed among clinical samples, there was not a 1 :1 
correlation between the expression of phosphorylated Smad2/3 
and NDRG2 in tissues. Thus, other factors may also be involved in 
TGF-p induced tumor suppressor gene expression and colon 
cancer progression, and some of these factors could be 
considered epigenetic alterations. 

Moreover, NDRG2 was found to be one of the intact TGF-p/ 
Smad signaling target genes induced in TGF-p-treated keratino- 
cytes. A cluster of sites located in the core region of the NDRG2 
promoter was further confirmed to contribute to inducibility by 
TGF-p in reporter assays. Although Spl was considered necessary 
for the basal transcription of multiple genes, we found that Spl 
binding at consensus sites in the core region of the NDRG2 



promoter was required for TGF-p-induced NDRG2 upregulation. 
Our data were similar to the observation of increased p15 lnk4B 
transcription in response to TGF-p, which is mediated through the 
interaction between Spl and the Smad2/3/4 complex. In 
addition, we also found that TGF-p induced a rapid reduction in 
c-Myc levels, thus abrogating the interaction between c-Myc and 
Miz-1, which was responsible for the transcriptional repression of 
NDRG2. Interference of the c-Myc/Miz-1 interaction with a 
dominant-negative construct, Miz-1 A(637-803), enhanced the 
TGF-p activating effect on the NDRG2 promoter. This finding 
further supports the hypothesis that TGF-p relieves the repression 
of c-Myc and increases Miz-1 transcriptional activity in NDRG2 
transcription. On the one hand, TGF-p/Smad signaling directly 
transactivated NDRG2, which is dependent on the recruitment of 
Spl to the NDRG2 promoter. On the other hand, as a second input 
for NDRG2 activation, TGF-p reduced c-Myc expression and 
abrogated the repression of the c-Myc/Miz-1 complex to allow 
for NDRG2 transcription (Figure 71; Supplementary Figure S9). A 
related mechanism was reported for p21 Clp1 as a target of TGF-p 
and c-Myc. 39 
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Figure 6. NDRG2 knockdown promotes TGF-p-induced EMT. (a) Phase-contrast photomicrographs of control cells and cells treated with TGF-p 
demonstrated that NDRG2 and p21 Qpl knockdown mimics the EMT phenotype in HT29 cells, (b) Circularity index of untreated and TGF-(31- 
treated cells, (c) Western blot analysis of E-cadherin and vimentin expression in total lysates from untreated cells and cells treated with 10 ng/ 
ml TGF-J3 for 48 h. (d, e) Migratory and invasive behavior of untreated cells and cells treated with 10 ng/ml TGF-p for 48 h. The control values 
were normalized to 1, and the data are expressed as the fold-change in treated cells, (f) The expression level of mRNAs encoding E-cadherin, 
N-cadherin, fibronectin, Twist, Snail and vimentin in control cells and cells treated with 10 ng/ml TGF-p for 48 h. (d-f) The data are presented 
as the means + s.d. from triplicate analysis. 



It is well established that most CRCs are resistant to the tumor 
suppressor effects of TGF-p and prone to TGF-p-induced EMT, 37 
which is associated with features of advanced disease including 
metastasis, resistance to chemotherapy and the generation of 
cancer cells with stem cell-like characteristics. 40 Specifically, the 
loss or abnormal intracellular localization of E-cadherin in colon 
cancer has been associated with advanced histological grade, 
metastasis and a reduction in disease-free and overall survival. 41 
Herein, we found that the induction of NDRC2 by TGF-p was 
completely abolished in CRC and that this induction accounted for 
the loss in anti-proliferative response. Furthermore, NDRG2 
overexpression remarkably attenuated TGF-p-induced EMT and 
decreased the migratory and invasive potential of HT29 cells in 
response to TGF-p. The reduction was stronger than that caused 
by p21 Clpl , one of the classic TGF-p responsive targets. Moreover, 
NDRG2 knockdown dramatically promoted EMT in HT29 cells 
mainly through the upregulation of key factors involved in TGF-p- 
induced EMT, including fibronectin, Twist, Snail and N-cadherin. It 
will be interesting to investigate how NDRG2 regulates these key 
EMT factors. Recently obtained evidence has demonstrated that 
NDRG2 significantly suppressed tumor metastasis by attenuating 
active autocrine TGF-p production in breast cancers, which 
suggests the need for an NDRG2-mediated anti-metastatic 
effects study. 42 There are data showing that NDRG2 was 
negatively regulated by TGF-p during the progression of 
hepatocellular carcinomas. 43 This observation may be due to 
impairment in the TGF-p/Smad signaling pathway or the 
activation of non-Smad signaling cascades (PI3K/AKT, p38MAPK 
and so on) in these cell lines in response to TGF-p. In addition, 
related evidence has shown that the enhancement of GSK-3P 
activity by NDRG2 overexpression causes proteasomal 



degradation of the Snail transcription factor and subsequent 
transcriptional regulation of EMT-related genes, 44 which is 
consistent with this study. Recent observations suggest that 
crosstalk occurs between the cell-cycle regulation and EMT 
pathways, 45 and p21Cip1 has been shown to attenuate EMT in 
Ras- and c-Myc-dependent breast cancer 46 We observed that p21 
loss is necessary for EMT promotion, whereas p21 overexpression 
has little effect on EMT, indicating that p21 is necessary for 
maintaining the epithelial phenotype but that it is not the sole 
reason. Despite this finding, the underlying mechanism requires 
further investigation. Nonetheless, our current data suggest that 
the NDRG2 tumor suppressor could inhibit TGF-p-induced EMT as 
well as cell invasion and migration in late-stage CRCs. 

To further understand whether NDRG2 controls the switch from 
TGF-p-mediated anti-proliferative gene responses to oncogenic 
potential in CRCs, we used clinical tumor samples to analyze the 
methylation status of the NDRG2 promoter, which is an important 
epigenetic regulation mechanism. Most colorectal tumors dis- 
played aberrant hypermethylation in the NDRG2 promoter 
compared with the hypomethylation in normal counterparts. 
Furthermore, in addition to being prone to carcinogenesis, NDRG2 
methylation also tended to occur in late-stage subjects, which was 
suggestive of the abrogation of TGF-p-mediated growth inhibition 
and increasing instances of oncogenic TGF-p properties in late- 
stage CRC. Thus, we hypothesized that repression of NDRG2 
resulted from promoter methylation, which occurs in late-stage 
colon cancer patients, and that this repression tipped the balance 
of the TGF-p role toward oncogenesis in CRC. 

To support this conclusion, the blocking intensity of promoter 
methylation was detected in A/D/?G2-hypermethylated colon 
cancer cells that overexpressed Spl. Strikingly, the induction 
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Figure 7. Promoter hypermethylation reduces NDRG2 expression and hinders the progression and metastasis of colorectal cancer, (a) The RNA 
and protein levels of NDRG2 in multiple colorectal cancer cell lines and the normal human intestinal epithelial cells (HIEC) were detected by 
comparative RT-PCR and western blotting, (b) Bisulfite sequencing analysis of the NDRG2 promoter in a panel of cell lines. (Top) Schematic of 
the NDRG2 promoter. The positions of the CpG dinucleotides are shown to scale by vertical lines. (Bottom) Each circle represents a CpG 
dinucleotide: open (white) circles denote unmethylated CpG sites and filled (black) circles indicate methylated CpG sites. Each row represents 
a single clone, (c) The methylation status of the NDRG2 promoter was analyzed by methylation-specific PCR in each of the 24 normal colon 
and colorectal cancer tissues. The white and black colors represent hypomethylation and hypermethylation, respectively. Each experiment 
was independently repeated three times, (d) The binding activity of Spl to methylated and unmethylated recognition sites in vitro. (Top) 
Double-stranded oligonucleotides were derived from the Spl recognition sites in the NDRG2 promoter. The methylated C (C m ) is shown in red. 
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Spl expression vector. Twenty-four hours post transfection, the cells were treated with 5-aza-dC at the indicated concentrations for 48 h 
before harvesting for western blotting and quantitative real-time PCR analysis. The data are the means ± s.d. from triplicate analyses, (f ) The 
immunoreactivity score (IRS) distribution of NDRG2 IHC staining among tumors of different invasion statuses, (g) The IRS distribution of NDRG2 
IHC staining between tumors of different lymph-node metastasis statuses, (h) The IRS distribution of NDRG2 IHC staining between tumors of 
different distant metastasis statuses, (i) Model of NDRG2 tipping the balance of oncogenic TGF-p via EMT inhibition in colorectal cancer. 



of NDRG2 by Spl was abolished in these cells primarily due 
to failed binding of Spl to the hypermethylated NDRG2 
promoter, whereas treatment with DNA methyltransferase 
inhibitor significantly restored this regulation. Therefore, the 
expression of TGF-(3-responsive tumor suppressor genes and 
their promoter methylation statuses should be valuable for the 
prediction of the role of TGF-p-mediated oncogenesis in CRCs. 



Additional studies on a much larger scale are needed to obtain 
more definitive conclusions. 

In summary, our data provide the first evidence that the NDRG2 
tumor suppressor gene is a new TGF-p responsive target and 
inhibits TGF-p-induced EMT in CRCs. The loss of NDRG2 expres- 
sion, which is attributed to promoter methylation, tipped the 
balance of the TGF-p role toward oncogenesis in late-stage CRC. 
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Table 3. Statistical results of IHC assay for NDRG2 in cohort of 260 
CRC 






NDRG2 
Negative Positive 




Total 


260 


1 84 


76 




Gender 








0.863* 


Male 


139 


99 


40 




Female 


121 


85 


36 




Age at diagnosis, years 








0.493" 


<60 


142 


103 


39 




>60 


118 


81 


37 




Depth of invasion 








0.002 b 


T1 


26 


14 


12 




T2 


62 


36 


26 




T3 


138 


104 


34 




T4 


34 


30 


4 




Lymph-node metastasis 








<0.001* 


Absent (NO) 


163 


98 


65 




Present (N1-3) 


97 


86 


11 




Distant metastasis 








0.006" 


Absent (MO) 


222 


150 


72 




Present (M1) 


38 


34 


4 




Abbreviations: CRC, colorectal cancer; NDRG2, N-Myc downstream- 
regulated gene 2; IHC, immunohistochemical. a P-value when expression 
levels were compared using Fisher's exact test. b P-value when expression 
levels were compared using Kruskal Wallis. 



Our observations elucidate how selective alterations in the TGF-p 
signaling pathway contribute to CRC development and progres- 
sion. These data will likely permit the identification of an 
additional subset of inherited CRC cases and provide novel 
opportunities for therapeutic intervention. 

MATERIALS AND METHODS 

Tissue samples and study cohort 

This study was approved by the ethics committee of the Fourth Military 
Medical University. All 50 patients who provided CRC and adjacent normal 
colon mucosa colectomy specimens and the 260 subjects in the CRC 
sample study cohort provided full consent for the study. The recruited 
patients were consecutively diagnosed with CRC in the Department of 
Gastrointestinal Surgery, Xijing Hospital, Fourth Military Medical University. 
Patients who met any of the following criteria were excluded: (1) received 
treatment before surgery that included neoadjuvant chemotherapy; 
(2) had insufficient specimens for tissue microarray preparation; (3) were 
diagnosed with gastrointestinal stromal tumor or lymphoma; (4) were 
diagnosed with additional cancers; or (5) refused to provide authorized 
consent. All of these clinical cancerous specimens were collected through 
surgery or endoscopy. All of the specimens were histologically diagnosed 
by the Department of Pathology, Xijing Hospital, Fourth Military Medical 
University. The clinicopathologic information for the remaining patients 
was entered into a database. For the study cohorts, study physicians who 
reviewed the CRC records and recorded the data into a database were 
blinded to the data. 

IHC detection 

IHC was performed on all CRC samples using the avidin-biotin-peroxidase 
method. All sections were deparaffinized in xylene and dehydrated using a 
concentration gradient of alcohol before endogenous peroxidase activity 
was blocked using 0.5% H 2 0 2 in methanol for 10min. After non-specific 
binding was blocked, the slides were incubated with a NDRG2 polyclonal 
antibody (1:200; ABNOVA, Jhongli, Taiwan) or p-Smad2/3 (1:100, sc-1 1769; 
Santa Cruz Biotechnology, Inc, Santa Cruz, CA, USA) in phosphate-buffered 



saline at 4°C overnight in a humidified container. Biotinylated goat anti- 
rabbit IgG (1:400; Sigma, St. Louis, MO, USA) was incubated with the 
sections for 1 h at room temperature and detected with a streptavidin- 
peroxidase complex. The brown color indicative of peroxidase activity was 
developed by incubating with 0.1% 3,3-diaminobenzidine (Sigma) in 
phosphate-buffered saline with 0.05% H 2 0 2 for 5 min at room temperature. 
The appropriate positive and negative controls were included in each IHC 
run. 



Evaluation of IHC staining 

Multiple tissue arrays were independently scored by two pathologists 
blinded to the clinicopathologic information and outcome of the patients. 
Staining was evaluated by scanning the entire tissue specimen under low 
magnification ( x 40) and then confirmed at higher magnification ( x 200 
and x 400). An immunoreactivity score system based on the proportion 
and intensity of positively stained cancer cells was used. The extensional 
standard used (1) the number of cells with positive staining (^5%: 0; 6- 
25%: 1; 26-50%: 2; 51-75%: 3; and >75%: 4) and (2) the staining intensity 
(colorless: 0; pallide-flavens: 1; yellow: 2; brown: 3). We multiplied the 
numbers scored in 1 and 2 from the extensional standard, and the staining 
grade was stratified as absent (0 score), weak (1-4 score), moderate (5-8 
score) or strong (9-12 score). All slides were interpreted by another 
pathologist unaware of the data. Specimens were rescored if the difference 
in score from the two pathologists was >3. Tumors with weak, moderate 
or strong immunostaining intensity were classified as staining positive ( + ), 
whereas tumors with no immunostaining were classified as staining 
negative ( — ). 



Cell culture and reagents 

The HaCaT human keratinocyte cell line, the Caco2 human colon cancer 
cell line and the HEK293 human embryo kidney cell line (ATCC) were 
maintained at 37 °C under 5% C0 2 in Dulbecco's minimal essential 
medium (Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal 
bovine serum and 1% penicillin/streptomycin. The human colon cancer 
cell lines HT29 and HCT116 (ATCC) were maintained in McCoy's 5A 
medium supplemented with 10% fetal calf serum. TGF-fSI (100-21) was 
purchased from Pepro Tech (Rocky Hill, NJ, USA). Cells were treated with 
10 ng/ml TGF-J31 for 20, 24 or 48 h, washed with phosphate-buffered saline, 
and harvested for downstream assays. 



MTT assay 

For cell growth and viability assays, cells were seeded into each well of a 
96-well plate with five replicates for each group at each time point. After 
the stated incubation time, 20 ul of MTT solution (5mg/ml) was added to 
each well, and after 4h of incubation, the medium was aspirated as much 
as possible without disturbing the formazan crystals. Then, 1 50 ul of DMSO 
was added to each well, and the plates were placed on a plate shaker for 
10 min. The OD values were determined at 570 nm using a Sunrise 
microplate reader (Tecan, Groedig, Austria). 



Real-time PCR 

For RT-PCR, cDNA was synthesized from RNA using the PrimeScript 
RT-PCR Kit (TAKARA, Dalian, China). We then subjected the cDNA 
(1-100ng) to PCR amplification. Quantitative real-time PCRs were 
performed using specific primers (Supplementary Table 1). The mean Ct 
values for the target genes were normalized to the mean Ct value for the 
endogenous control GAPDH. The ratio of mRNA expression of target gene 
versus GAPDH was defined as 2 (-ACt). 



Plasmids 

Human NDRG2 promoter truncated regions were previously cloned into 
the pGL3-basic vector (Promega, Madison, Wl, USA) to generate NDRG2I 
luciferase reporter plasmids. 20 The NDRG2 promoter mutants were 
generated using the mutation primers listed in Supplementary Table 1. 
pcDNA3.1 ( + )-c-Myc was generously provided by Dr Michael Cole (Lewis 
Thomas Laboratory, Princeton University, Princeton, NJ, USA). pcDNA3.1 
( + )-Miz-1 was kindly provided by Dr Ben-Zion Levi (Department of 
Biotechnology and Food Engineering, Technion, Haifa, Israel). 
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Western blotting 

Cells were collected and lysed in lysis buffer. Protein concentrations were 
measured using the BCA protein assay kit from Thermo Scientific 
(Waltham, MA, USA). Forty micrograms of protein were separated by 
SDS-PAGE and transferred onto Hybond ECL nitrocellulose membranes. 
Western blot analyses were performed using the following primary 
antibodies: anti-WDRG2 mouse monoclonal antibody (H00057447-M03), 
which was purchased from ABNOVA; anti-Sp1 rabbit polyclonal antibody 
(sc-59), anti-E-cadherin polyclonal antibody (sc-7870) and anti-vimentin 
mouse monoclonal antibody (sc-373717), which were purchased from 
Santa Cruz Biotechnology, Inc.; anti-p21 Clp1 rabbit polyclonal antibody 
(#2947), which was purchased from Cell Signaling Technology (Beverly, 
MA, USA); and mouse monoclonal anti-human fS-actin (1:5000 dilution), 
which was purchased from Sigma. Secondary antibodies were chosen 
according to the species of origin of the primary antibodies and detected 
using enhanced chemiluminescence (Pierce, Rockford, IL, USA) or the 
Odyssey Imaging System (Li-Cor Biosciences, Lincoln, NE, USA). 



Bisulfite-PCR and sequencing and methylation-specific PCR 
Total genomic DNA was isolated from multiple cell lines and the tumor and 
normal colon epithelium samples using a DNA extraction kit according to 
the manufacturer's instructions (Tiangen, Beijing, China). Two micrograms 
of genomic DNA was subjected to a denaturation/incubation cycle in 
sodium bisulfite solution, as previously described. 47 For bisulfite 
sequencing analysis, 10 clones were sequenced for each cell line. For 
methylation-specific PCR, the intensity of the methylated and loading 
control bands was compared with the base pairs of the amplicons, and the 
cycles of individual PCRs were considered. When the relative ratio of 
methylated DNA to the total amount of DNA was greater than 50%, the 
methylation status was considered to be high, whereas it was considered 
low when the ratio was less than 50%. Methylation was considered absent 
in samples with no obvious bands (see Supplementary Table 1 for primer 
sequences). 



Luciferase reporter assays 

Cells were transfected at ~ 70-80% confluence using Lipofectamine 2000 
(Invitrogen) according to the manufacturer's instructions. Twenty-four 
hours after transfection, the cells were treated with 10ng/ml TGF-p for 
20 h. Luciferase assays were performed using the dual luciferase reporter 
assay system (Promega). The values are expressed as the means ± s.d. of at 
least three independent experiments. 



Chromatin immunoprecipitation 

ChIP analysis was performed using the ChIP Assay kit (Upstate 
Biotechnology, Charlottesville, VA, USA). The cells were crosslinked using 
1% formaldehyde for 10min at 37 "C and then washed, lysed and 
sonicated to generate 200-500 bp DNA fragments. The samples were 
precleared using 60 ul of salmon sperm DNA-protein A-agarose and 
subsequently incubated at 4 °C overnight with 2 ug of HA or Myc antibody; 
rabbit IgG was used as a control. Immunocomplexes were recovered, 
washed thoroughly and eluted with ChIP elution buffer. Following the 
reversal of crosslinks at 65 °C for 4h, the DNA from each sample was 
extracted using phenol/chloroform, precipitated with ethanol and then 
used as a template for PCR amplification. PCR primers were used to detect 
the NDRG2 promoter (Supplementary Table 1). 



In vitro migration and invasion assays 

The migration and invasion of CRC cells were examined using 
polycarbonate transwell filters containing 8 urn pores (Becton Dickinson 
Labware, Franklin Lakes, NJ, USA). After treatment, the cells were seeded in 
serum-free media on the upper side of a transwell chamber that was either 
uncoated, for the migration assay, or coated with Matrigel (BD Biosciences, 
Bedford, MA, USA), for the invasion assay. The cells were allowed to 
migrate toward media containing 10% fetal bovine serum for 24 h. After 
the incubation period, the cells on the lower side of the membrane were 
fixed, stained with crystal violet, and counted. The migration and invasion 
indices were calculated as the mean number of cells in 1 0 random fields at 
x 20 magnification. 
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Circularity index analysis 

Circularity index analysis of control cells and cells treated with 5 ng/ml of 
TGF-p1 for 48 h was performed. For each condition, the circularity of > 100 
cells in at least two separate areas was determined, and their average ± s.d. 
was determined. *P<0.05 compared with control cells. 

Statistical analyses 

The data are expressed as the means ± s.d. Statistical analyses using 
Student's T-test for independent groups was performed using the SPSS 
16.0 software package (SPSS Inc, Chicago, IL, USA) for Windows. 
Associations between NDRG2 expression and categorical variables were 
analyzed by the Mann-Whitney Otest or the Kruskal-Wallis test, as 
appropriate. *P<0.05 was considered as statistically significant. 

ABBREVIATIONS 

TGF-p, transforming growth factor |3; NDRC2, N-Myc downstream- 
regulated gene 2; EMT, epithelial-mesenchymal transition; 5-aza- 
dC, 5-aza-2'-deoxycytidine 
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